Non equilibrium electronic distribution in single electron devices 
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The electronic distribution in devices with sufficiently small dimensions may not be in thermal 
equilibrium with their surroundings. Systems where the occupancies of electronic states are solely 
determined by tunneling processes are analyzed. It is shown that the effective temperature of the 
device may be higher, or lower, than that of its environment, depending on the applied voltage 
and the energy dependence of the tunneling rates. The I-V characteristics become asymmetric. 
Comparison with recent experiments is made. 
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In small devices, the coupling between electrons and the lattice is suppressed. At sufficiently low temperatures, 
only long wavelength accoustical phonons are excited. If the electrons are localized in a region much smaller than the 
wavelength of the phonons, the interaction between the electrons and the phonons is significantly reduced As a 
result, the electron temperature may be different from that of the surrounding medium. 

Usually, the electron temperature in small devices tend to be higher than that of the environment, because of 
dissipation at the device (see also ||^), via shake-up processes. This effect can explain a number of experiments 
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In the following, we analyze the electronic distribution in a small device, when it is controlled by the electronic 
exchange with the external leads. The general equations which describe the single level occupancies were first discussed 
in |^,|| . These authors considered mostly the inluence of a non negligible single level spacing in the I-V characteristics 
of semiconductor nanostructures. In these systems, the effects of the energy dependence of the transmission through 
the barrier, or of the density of states in the external leads is not important. We will generalize the previous work 
to situations where energy dependent tunneling rates also contribute to modify the electronic distribution. As it will 
be seen below, an energy dependent tunneling rate, T(e), can modify significantly the electronic distribution in the 
central electrode of a single electron device, if T'(e)/T(e) is comparable to the inverse of the temperature at which the 
device is operated. Such a situation can be realized when an electrode has a strongly dependent density of states, or if 
the tunneling process takes place close to the top of a barrier. The first case is relevant to the experiments reported in 
0, where one of the electrodes is made of graphite, and to the experiments presented in where the electrodes are 
made of superconducting Al, and tunneling processes take place near the gap edges. Asymmetric I-V characteristics 
have also been reported in . Some features of these experiments are also consistent with the work reported here. 

A particularly interesting case is presented in The observed Coulomb staircase can only be fitted by the 
orthodox theory |ll|-|l3||, if an effective temperature of ~ 170K is assumed, although the experiment is performed 
at room temperature. A large temperature difference between the central electrode in the device used in and its 
surroundings cannot be ruled out. The central island is isolated by a Langmuir-Blodgett film from the rest of the 
system. A simple estimate of the temperature difference can be made by assuming that the energy dissipated in the 
circuit has to be carried away through the film. Then: 
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where k is the thermal conductivity of the film, L is the linear size of the electrode, T^nv is the temperature of the 
environment, and Teiec is the temperature of the electrode. / and V are the intensity and the voltage acting on the 
device. Typical values of k for organic films, such as Mylar, are ~ IQ-^ W / ( m K ). The size of the electrode is 
L 20A. Using the values for / and V reported in 0, we find that T^nv — T^iec ^ lO^K. Note that the temperature 
gradient supported by a good thermal insulator, such as an organic film, is much larger than the one expected in a 
system where heat can be dissipated through metallic regions (k~1W / (mK)). 

A different way to manipulate the electronic distribution in a mesoscopic device has been discussed in Our 
analysis differs from that in |^ because we will study systems in which the electronic distribution is not in thermal 
equilibrium. 

The orthodox theory of Coulomb blockade [pdj-p^ analyzes the transport processes in small capacitance devices by 
means of rate equations. These equations assume that the different elements of the device are in thermal equilibrium. 
The rate equations allow us to obtain the (time dependent) number of electrons in each individual electrode. We will 
relax the hypothesis of electronic thermal equilibrium, as in ||J^. We assume that transport processes determine, not 
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only the evolution of the charge, but the fluctuations in the occupancies of individual levels as well. The influence of 
other effects, such as electron-electron collisions or the coupling to the lattice, will be discussed later. 

We study a small island, coupled by two junctions to leads connected to a battery which maintain a voltage 
difference, V. Each junction is characterized by a capacitance and a resistance, Ci^2 and i?i,2- The electronic 
temperature in the leads, T, is fixed by the cryostat. We will assume that the charge in the island fluctuates between 
two states, with N and + 1 electrons. The island has an effective charging energy of e^/(Ci + C2) (l^. When the 
island has charge iV, the voltages at each junctions are [T^ : 



Ci + C2 Ci + C2 

F2 = -^^ + y^^ (2) 
^ C1 + C2 C1 + C2 ^ ' 

These voltages determine the tunneling rate of electrons to and from the leads. 

We label the probability that the island is in a state with N electrons in the levels 11,12, ■■■,iN as ni-^^i^^,,,^ij^ , and 
we use a similar convention for the cases with N + 1 electrons. The probability for an electron from lead 1 to jump 
into a state of energy e in the island is proportional to the occupancy of a given state at that energy, n(Vi + e), 
times the transition rate across junction 1, ti{Vi + e), where n(e) is the Fermi-Dirac distribution at temperature T. 
For simplicity, we include density of states effects in the definition of the tiS,. Analogously, the probability that an 
electron of energy e leaves the island through junction 1 is proportional to f — n{Vi + e) times ti{Vi + e). We define, 
in this way, the quantities: 

pT{e) = h{Vi+e)n{Vi+e) 
pT\e) = h{Vi+e)[l-n{Vi+e)] 
p':^{e) = t2{V2 + e)n{V2 + e) 

pT\e) = t2{V2 + e)[l-n{V2+e)] (3) 
Then, the probabilities of finding a given electronic configuration in the island obey: 
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lfc=Jl...ljV + l 
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These equations are a straightforward generalization to the ones formulated in ||^,^ to the case of energy dependent 
tunneling rates. They admit the stationary solution: 



where C is a normalization constant and: 

■'^ ' pr*(e)+pr*(e) 

Charging effects are built in into eq.(|5|) through the dependence of the p^'s on Vi and V2. 

Each factor in the product in eq.dS) can be interpreted as a Boltzmann weight which determines the occupancy 
of the corresponding level. The previous scheme can be extended to situations where more than two charge states 
are involved. Then, each pair of charge states, N and N + 1, determine a set of equations of the form (Q) [^. The 
quantities and depend on N and N + 1. A simple solution, like (^), is no longer possible. A closed expression 
can be obtained in terms of sums over all charging processes which lead to a charge state N, with occupancies {ik}, 
from the minimum charge state included in the model, Nmin- Finally, in the absence of charging processes, the 
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voltages Vi and V2 are independent of the charge state of the central islands. The solution of the rate equations (Q) 
is given by (H). The calculation of Vi and V2, however, becomes more complicated, as they are determined by the 
junction resistances, and a self consistency loop is required. 

In the following, we will describe the electronic distribution by means of /(e), (|^). This scheme is valid for voltages 
in the flat portions of a Coulomb staircase, when only two charge states are involved. Close to the steps, the island 
fluctuates between three charge states, N,N + 1 and N + 2. The chemical potential for state + 2 is slightly above 
the chemical potential of one of the leads. Then, the rates and for processes which take the island in or out 
of state + 2 depend more weakly on energy than the processes which relate states N and + 1 . In this case, a 
description of the occupancies in the island by a function /(e) which depends only on the rates between N and 1 
is adequate. 

The distribution described by /(e) is not in thermal equilibrium. This is reflected in the fact that the quantity 
defined as /3(e) = —d\og[p{e)]/de is not independent of e, although the variations of /3(e) will be less pronounced than 
those of p{e). If /3'(e) <^ /3^(e), an approximate thermal distribution can be defined in an energy range ~ /3~^ around 
e. The effective temperature is (3^^. In the central electrode, unless /3 has a very unusual dependence on e, te lowest 
lying states will be fully occuppied, and the highly excited ones, empty. Hence, if (3' <SC /3^ around the equilibrium 
Fermi level, the new distribution is completely specified by 1/ {kBT^f f) = (3{ep). T^ff, defined in this way, gives the 
best fit to the level occupancies within the island by a thermal distribution. In principle, it depends on e. 

The following properties of Tef / are easy to prove: 

(i) T,fj = T if y = 0, for any t,^2{e). 

(ii) Teff —^Ti^tl and t2 are independent of energy, and ^2 ^ ti, or ti ^ ^2- The thermalization of the central 
electrode in a very asymmetric junction was discussed in [ p^ . 

(iii) Tf^ff = T, if ti = ^2- The lack o corrections to the equilibrium distribution in a symmetric junction, to lowest 
order in the applied voltage, was discussed in . 

(iv) li V ^ Q and ti , t2 are independent of energy, T < T^ff. T^ff reaches a maximum for Vi < e < ¥2- If 
Teff — !■ 00, all electronic configurations with a given TV are equally likely. This situation may have been observed 
experimentally jll. Note, however, that, if Te,, -. 00 , /?' ^ /3^ . The electronic distribution cannot be well described 
by Teff in an energy range comparable to Te// itself, in this case. 

We now consider the situation when the ti's depend on e. We study the case t2 ^ ti, which is a neccessary condition 
for the observation of a Coulomb staircase Then, the leading effects associated to t[,t2 give: 
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n{V2 + e) _ 1 - n{V2 + e) 
n{Vi + e) 1 - n{Vi + e) 



— < min[n(Vi + e), 1 - n{Vi + e)] 
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ksTeff tl tl 

^ - l-n(l/i+e)«^ (7) 
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The most favorable case for cooling occurs when 7i(e + Vi) <C ^2/^1- Then, T^ff ^ nim[T,t[/{kBti)]. If the junctions 

are tunnel barriers, ti{ei) oc r , where C is a constant of order unity, Vq is the height of the barrier, m is 

the mass of the electron and I is the length of the barrier. Then, 
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Typical values for m, Vb and I give t'^/ti ~ 1 (eV)~^. Density of states effects, in a normal metal, give modulations 
to the tiS of similar magnitude. Hence, the e dependence of the tunneling rates should not induce large effects in 
conventional SET's, unless the Fermi surface of the central island lies close to the top of one of the barriers. 

Significant changes in the T^ff^s, are expected if the tunneling rates, or the density of states in the substrate, 
depend strongly on energy. This is the case in the experiments reported in The Fermi energy of the central 
island is aligned with the edge of the gap of one of the superconducting electrodes. Thus, the process studied here 
can contribute to the cooling observed in |^ . A different situation is the setup used in . One of the substrates is 
graphite, which is a semimetal. The density of states per atom^near the Fermi level, increases roughly linearly with 
energy, ^(e) - |e - A|/e§, where eo - 4eV, and A ~ IQ-'^eV Hence, keT^ff ~ D / D' ~ \Vi - A|-i. Near the 
steps of the Coulomb staircase, the Fermi level of the island is close to that of the substrate, so that Vi measures 
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the distance to the step edge. The effective temperature can be well below that of the cryostat. Note that, for a 
tunneling rate which depends linearly on energy, /3'(e) ~ so that the characterization of the distribution by a 

single parameter, T^ff, is only an approximation. The overall effect on the Coulomb staircase, however, is a significant 
sharpening of the steps, in agreement with 

The general situation of energy dependent barriers, or density of states, is schematically shown in fig.(0). Intuitively, 
if the value of ii(e) rises sufficiently fast with e, and electrons flow into lead 1, the high energy electrons are blown away 
from the island. If the applied voltage is reversed, hot electrons are pushed into the island, where they are blocked by 
the high resistance junction, 2. The I-V characteristics will be highly asymmetrical. A population inversion, in which 
the high energy states show larger occupancies than the low energy ones is also possible. This situation corresponds 
toTe// <0. 

So far, we have assumed that the electron distribution within the island is entirely determined by the tunneling 
processes to and from the leads. At the beginning of the article we presented an estimate of the efects o heat conduction 
across the surrounding medium on the temperature of the central electrode. We need also to consider electron-electron 
interactions, among the electrons in the central electrode and those in the leads, and within the electrode itself. In 
the orthodox theory of Coulomb blockade, tunneling processes are assumed to involve no dissipation. The relaxation 
of this approximation leads to shake up effects, and heating Q]. These processes impose a lower bound on the island 
electronic temperature, T^u, which is a fraction of the Coulomb energy, Ec. Taking fcsT,,, - E^/S § and E^ ~ 0.12eV 
§, we find Tsu ~ 170K. 

Finally, electron-electron interactions within the island will tend to thermalize the electronic distribution. The 
thermal energy in the island is proportional to T^ff, taken at the Fermi level. Hence, the equilibrium temperature 
will not be very different from this value. 

One of us (N. G.) is thankful to M. Jonson and R. I. Shekter for helpful discussions. 
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FIG. 1. (a) Combination of external voltage and barriers for electron cooling (high energy electrons leave the central island 
faster than the low energy ones), (b) Combination of external voltage and barriers for electron heating (high energy electrons 
arrive to the island faster than the low energy ones). 
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